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Like the even chirality correlation functions of, say, the two vector currents, one can also consider 
odd chirality correlation functions to write thermal QCD sum rules. They contain fewer non- 
perturbative corrections, at least to the leading order. Here we write such a sum rule for the 
correlation function of vector and tensor 'currents'. The odd and even chirality sum rules are taken 
together to evaluate the effective parameters of the p meson to second order in temperature. To 
within errors, the results are consistent and reproduce the absence of shift in the p meson mass to 
this order. 



I. INTRODUCTION 

The vacuum QCD sum rules jjj, when extended to finite temperature ||, provide a simple means to study the 
properties of the QCD medium. Below the critical temperature, they relate the temperature dependence of hadron 
parameters to the thermal average of local operators (|-[| . 

Generally the sum rules are written for correlation functions of two currents, which have the same chirality, eg, of 
two vector or two axial vector currents. But one may also consider correlation functions of two quark bilinears of 
opposite chiralities ||. Then their operator product expansions will be dominated, in each dimension, by the operators 
of odd chirality, the contributions of even chirality ones being highly suppressed by a factor of (small) quark mass 
in their coefficient functions. As a result the leading nonperturbative corrections will be given by a single operator, 
namely, the quark condensate; new odd chirality operators appear only as non-leading terms. 

Here we write down the sum rule for the correlation function of the vector current and the tensor 'current'. Then 
both the present vector-tensor sum rule and the two vector-vector sum rules derived earlier get contributions from the 
same set of intermediate states with J PG = 1 h . We consider the subtracted sum rules, obtained by subtracting out 
the vacuum sum rules from the corresponding thermal ones. The only intermediate states with significant contributions 
are then the non-resonant 2tt state and the /9-resonance. The three sum rules form a convenient set for evaluation . 

We confine our evaluation at low temperature T to leading order, which is T 2 in the chiral limit. Generally the 
thermal sum rules are complicated by the presence of Lorentz non-scalar operators, in addition to the Lorentz scalars 
contributing to the vacuum sum rules j8j. But to order T 2 non-scalar operators cannot contribute to these sum rules. 
Thus to this order, the contributing operators are the same as at zero temperature, their vacuum expectation values 
being replaced by the thermal averages. 

In Sec II we describe the kinematics and present the result of operator product expansion. In Sec III we saturate 
the sum rules with 2ir and p states. In Sec IV the sum rules are evaluated for the parameters of the p meson to order 
T 2 in the chiral limit. Finally in Sec V we discuss the different aspects of the thermal sum rules and compare our 
evaluation with earlier ones. 



II. SUM RULES 

We restrict here to the better known non-strange channels of unit isospin. The quark bilinears in this channel are, 

S a (x), V-(x),T^(x), A^(x), P a (x) = q(x)(l, lfl , a uv , ltl75>l5 ) ^-q(x), 

q(x) being the field of the u and d quark doublet and r a the Pauli matrices. Note that S a (x), P a {x) and T^(x) are 
odd and V^(x) and A a ^ (x) are even under 75 transformation on the quark field. 

As an important example of odd chirality correlation functions, consider the thermal average of the time ordered 
(T) product of the vector current and the tensor 'current', 
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T^{q) = i J d^xe* * (T (F;(x)T^(o))) . (2.1) 

Here the thermal average of an operator O is denoted by (O), 

(O) = Tre-^O/Tre- 011 , 

where H is the QCD Hamiltonian, /3 is the inverse of the temperature T and Tr denotes the trace over any complete 
set of states. 

As usual, it is convenient for kinematics to restore Lorentz invariance by introducing the four- velocity it^ of matter. 
Then we can build the Lorentz scalars, u) = u • oand q = y^ui 2 — q 2 , representing the time and the space components 
of q^ in the matter rest frame, it M = (1,0,0,0) |9|. We now choose the three independent kinematic covariants as, 

-f/ict/3 = Vnaqp — Vvpqu, 

Qfiap = q^{q a up - qpu a ) - q 2 (r)!_ ia up - r]^u a ), 

R^af3 = u^qaUp - qpu a ) ~ U^fiallp - r]^u a ). (2.2) 

The kinematic decomposition now reads, 

VaM = ^ ab (P^T! + Q mP T 2 + R^n), (2.3) 

where the invariant amplitudes I! 2,3 are functions of w and q 2 . In all computations, however, we shall revert back to 
the matter rest frame. 

Only the amplitude T\ survives at zero temperature. As we shall see below, at finite temperature the leading 
contributions (to order T 2 ) are also contained in T\. So, working to the leading order, we have only to write the sum 
rule for this amplitude. 

The advantage with odd chirality correlation functions becomes evident from an enumeration of local operators. 
The unit operator corresponds to the perturbative result. The non-perturbative power corrections begin with op- 
erators of dimension three and four. At dimension three, we have qq and q^,q belonging to odd and even chirality 
respectively. (The operator qi/iq actually cannot contribute for zero chemical potential.) At dimension four, we have 
trG^G^, u^G^u" and u^O^u". Here the gauge field strength = g^-G a ^, X a (a = 1, • • -8) are the SU{3) 
Gell-Mann matrices and g is the QCD coupling constant. The operators Qv§ are the energy momentum tensors 
for the quarks and gluons. (Note that in the matter rest frame, the latter two are just the energy densities. Also 
the operator qu ■ Dq where is the covariant derivative, is of odd chirality, but it can be reduced to rhqi^q by 
using the equation of motion for the quark field.) All these operators of dimension four are of even chirality. Thus 
up to dimension four, only the operator qq can contribute significantly to the power correction for an odd chirality 
correlation function. At dimension five, there is only one Lorentz scalar operator, qa^G^q. In addition there are 
several Lorentz non-scalar operators contributing at finite temperature |Tof . 



There is a simple configuration space approach 



of all operators, both scalars and non-scalars 12 
for the product under consideration, 



1 1 1 to the operator product expansion giving the Wilson coefficients 
12]. Using this method and restricting to scalar operators, we get 



tv«( X) tuo) - n M - ^p^ 1 + i j^w Uu + 21b (^) 0s + -\ (2 - 4) 

where O5 is the dimension five operator, O5 = ua^G^u, and the dots represent operators of still higher dimensions. 
We have assumed SU (2) flavour symmetry; m is the degenerate mass of u and d quarks and uu — dd — \qq. The 
Fourier transform gives for space-like momenta (Q 2 = —q 2 > 0), 

f 3m 11 1 

"> ^(V^qn - Vtfq«) {-^log(Q 2 /fi 2 ) - -^{uu) + — ^- (0 5 ) + • • • j . (2.5) 

Here u (~ 1 GeV ) is the renormalisation scale. 

We wish to include the renormalization effects on the operators and uu. The coefficient C(Q 2 / /j?, g(p)) of uu 
satisfies |l4| , 

3d \ 
^ + P{9)Q- g + 7i - 72 j C(Q 2 /f, 2 ,g(p)) = 0, (2.6) 
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where (3{g) = ~b j£pp , b = 9. The anomalous dimensions 71 and 72 of the operators T® v and tut can be easily 

2 2 

calculated to give 71 = d\ rfcs and 72 = d,2j£pj? with <ii = 8/3 and e?2 = — 8. The solution to (2.6) may be written as 

C(Q 2 /p 2 ,g(v)) = a(Q 2 )C(l,g(Q 2 )), 
where C(l,g(Q 2 )) is the lowest order result obtained above and 

^ 2 ,_ f log(Q 2 /A 2 ) \ (dl - d2)/2b 
a{Q) -\log{^/k 2 )) 

The strong interaction scale is A ~ 200 MeV. 

The other element needed for the sum rule is the spectral representation for the correlation function. The T product 
at finite temperature has the so-called Landau representation in the variable qo at fixed |<f| Jl5| ] . At points on the 
imaginary axis (q 2 , — — Qo>*3o > 0)j the representation for the invariant amplitudes are, up to subtractions, given by 



T la ,z{qlAq\) = [ d<7o^7 
Jo % 



P N lt2t3 (q' 0> \q\) 



9o +<2o 



(2.7) 



where 

iVl,2,3(<?0, \q)) = ^^mTi^s^o.kDtanh^o^). 



III. SATURATION 



In the channel under consideration, the p- meson dominates the absorptive part. To find this contribution we write 
the matrix elements of the currents, 

(0\V^\p b (k)) = 6 ab F p m p e p , (0\T^\p b (k)) = i8 ab G p {e a k p - tpk a ). (3.1) 

Here F p and G p are two coupling constants and m p and e a are the mass and the polarization vector of the p-meson. 
F p is measured by the electronic decay rate of p°, F p = 153 MeV [^6|. Though the value of G p is not available 
directly from experiment, it can be obtained from one of the U(6) symmetry relations for the wave functions of a 
quark-anti-quark pair JlTj . Defining the pion decay constant F n by 

<0|^|vr b (fc)) (3.2) 

this relation is 

G p = \{F p +F n ). (3.3) 

As Fn = 93 MeV, we get G p = 123 MeV. This value of G p is also obtained from the QCD sum rules for the correlation 
function of two tensor 'currents' ||. Comparing (3.1) with the matrix element of the p-meson field operator p^(x), 
(0\p p \p b ) = 5 ab e^, we get the operator relations, 

V«{x) = m p F pP *(x), T^(x) = G p (B a pl - dp ft). 

The absorptive parts are now given essentially by that of the p-meson propagator at finite temperature. Working 
in the real time formulation of the thermal field theory Jlq| , it is the 11-component of the 2 eg) 2 matrix propagator. 
We get 

N^q 2 ) = m p F p G p 6(q 2 - m 2 p ), N 2 (q 2 ) = N 3 (q 2 ) = 0. (3.4) 

This calculation must be interpreted as one in the effective field theory at finite temperature, where loop corrections 
make the parameters m p ,F p and G p temperature dependent. (At finite temperature, each of the particle-current 
coupling constants in (3.1-2) has, in general, different temperature dependence for the time and space components of 
the currents [l9| . But this bifurcation takes place only in orders higher than T 2 . So it does not concern us here.) 
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At lower energies there is the contribution of non-resonant two-pion state. Although small compared to the p-meson 
contribution, it describes the interaction of the currents with the pions in the heat bath and may assume importance in 
the difference sum rules we shall consider below. We find this absorptive part by writing the field theoretic expression 
for the pion loop at finite temperature. The pionic content J^(x) of the quark vector current V®(x) is given, to lowest 
order, by 

W - W = e abc cb b (x)d^ c (x). 

However the pionic content S^g(x) of the quark tensor 'current' T®o(x) is not immediately known, as it is not a 
symmetry current. From its index structure it may be written as 

T^{x) - S a a/3 (x) = ce abc d a c( > h (x)d^ c (x), 

where we determine the constant c by comparing the two divergences, 

d a T^(x) = 2rhV$(x) + ■■■, d a S^(x) = cm*J^(x) + • • • , 

the dots standing for higher derivative terms. Like V£(x) and J^(x), we can also identify T£ g (x) and S®p(x), giving 
c = 2m/m 2 = — F%/(0 \uu\0), on using the Gell-Mann, Oakes and Renner relation poj . 

With the pionic version of V?(x) and (x), it is simple to evaluate (2.1) to lowest order as, 



T$ 3 (q) = -c6 ab I t|4(2/.' ( /» ; ,(^ 1 /,';- 1 y i /-iA Ji( /,-)A. 1! .(/,- - H ). ,».:,) 



where An is the 11- component of the thermal pion propagator. Its absorptive part can be calculated in the same 
way as for the vector- vector correlation function As we need to write the sum rule for the amplitude T\, we quote 
the results for this amplitude only. In the time- like region (superscript +), 
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N i=J^I dx(v 2 (q 2 )-x 2 ){l + 2n((\q\x + q )/2)}, for q 2 > 4m*, (3.6) 



while in the space-like region (superscript — ), 

Ar = dx(v 2 (q 2 )-x 2 ){n((\q}x-q Q )/2)-n((\q]x + q )/2)}, for q 2 < (3.7) 

Here we have defined the functions v(z) = — 4m 2 j z and n(z) = {eP z — 1) — 1 . 
For|q| — > 0, they reduce to simple expressions. In the time-like region, 

N+ = C -^^{l + 2n( qo /2)}. (3.8) 

In the space-like region, q 2 = A|<f| 2 , with < A < 1. Thus in the limit |g| — > 0, A-f and its contribution to the spectral 
representation (2.7) are zero. 

We can now write a spectral sum rule by equating at a space-like point the operator product expansion (2.5) to 
the spectral representation (2.7) with p and 27r-contributions. Taking Borel transform as usual to get rid of any 
subtraction constant and to improve convergence of the spectral integral, we arrive at the thermal QCD sum rule, 

1^2 roc 

m p (T)F p (T)G p (T)e-<W M2 - f-_ / dsse-°l M \\s){l + 2n{^/2)) 

967^(0 | uu | 0) J 4m 2 

= |^M 2 - a{M 2 ){uu) + (0 6 ). (3.9) 



IV. EVALUATION 



As already stated, we shall evaluate sum rules to order T 2 in the chiral limit. Since the continuum integral is 0(T 4 ), 
the odd chirality sum rule (3.9), after subtraction of the corresponding vacuum sum rule, simplifies to 
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m p {T)F p {T)G p {T)e~<^ M2 - m p F p G p e- m > i2 = -a(M 2 )(uu) + ^(0 5 ). (4.1) 

Here the bar over the thermal average indicates subtraction of the corresponding vacuum expectation value. 

For the sake of consistent evaluation in a closed framework, we augment this sum rule with the two even chirality 
sum rules for the vector- vector correlation function derived earlier j^] . In the chiral limit and omitting terms of order 
higher than T 2 , they become 



4 M 4 



(4.2) 



Fl{T)ml{T)e-^ T ^ - F 2 p m 2 p e-<l^ = | ^ (4.3) 

where 

MM 2 ) = — ^ M^e^ 2 )^-^-^) (4.4) 

and Oq is the four-quark operator, 6 = ot s 97^75 A a r 3 q (77^75 X a r 3 q 21 1. This operator was ignored in 

At low temperature the heat bath consists primarily of dilute pion gas. The thermal trace can then be approximated 
by the vacuum and the one pion state. The pion matrix element of the operator can be evaluated by using PCAC 
and soft pion technique ||. One gets, 

T 2 \ ( T 2 \ ( T 2 

A. (Ok) =(0 ft Oi l . 



(uu) = (0 I uu I 0) II - ^ j , (0 5 ) =(0|O 5 |0) (1--^), (0 6 > = (0 I O e I 0) ^1 - ^ j . (4.5) 

where the vacuum expectation values are all known |^2| Q . 

(0 I uu I 0) = -(225MeV) 3 , (0 | 5 I 0) = m 2 o (0 \ uu | 0), (0 | O e | 0) = 6.5 x l(T 4 GeV 6 . (4.6) 
with m 2 = l.OGeV 2 . If we write 

m p (T) = m p (l + a^, F P (T) = F p [\ + b^j , G P {T) = G p (l + c^j , (4.7) 

the sum rules (4.1-3) predict the values of the constants a, b and c as 

b= -e 



18 4M 4 V m 2 p _ 

If m 2 \ K ? ( m 2 



K/^J_L h_^£ +— 1 + -^ H (4 9) 

1 18 V M 2 I 4Af 4 V M 2 ' 



m » „ / ,,,9, ml \ K 2 ( m 2 2M 2 



c 



where 



<o | »"|Q) = JL (o l o 6 1 o>. 

4m pJ F 2 (l + F 7r /F p )' 1 1 ' 



It should be noted that the parameters m p , F p and F v in eqn. (4.8-10) refer to their chiral limits and not the 
physical values. There does not appear to exist any calculation of F p to this limit in the literature p5| . Under this 
circumstance, let us evaluate a, b and c with the physical values of these parameters as quoted in Sec. Ill and hope 
that the results will not be affected much. 
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(GeV 2 ) 


a 


b 


c 


0.8 





-0.086 


-0.21 


1.0 


0.020 


-0.069 


-0.24 


1.2 


0.035 


-0.063 


-0..26 


2.0 


0.049 


-0.057 


-0.30 


4.0 


0.054 


-0.056 


-0.33 







Table 1: Coefficients of T 2 in m p (T), F p (T) and G P {T) at different values of M 2 . 



Table I shows the evaluation of the constants a, b and c for different values of the Borel parameter M 2 over the 
range 0.8 < M 2 < 4 in GeV 2 . It is seen that as M 2 increases, the values change but rather slowly. Observe that at 
M 2 = 0.8Gey 2 , one gets a = 0, reproducing the result that the hadron masses do not shift to order T 2 p6| , p7[ . This 
suggests that we take the central values of all the constants at this value of M 2 in the prediction of our sum rules, 

= 0, 6 =-.086, c=-.21, (4.11) 

with an error of about 25%, estimated from their variations in the range of M 2 considered. 

Eq(3.3 ) relating G p to F p and F„ is used above only at zero temperature. Being a symmetry relation one expects 
it to be valid also at finite temperature, 

G p (T) = ~(F p (T)+F n (T)). (4.12) 
This equation provides another determination of c, if we take F p (T) as determined above and 

i? ^ T ) =F -( i -iS^) ( 4 - ia ) 

as obtained from chiral perturbation theory. It gives c = —.14, in fair agreement with c = —.21 obtained above. 



V. DISCUSSION 



We derive here an odd chirality thermal QCD sum rule. More such sum rules can be written by considering other 
correlation functions of odd chirality, for example, of the axial vector and pseudo-scalar currents. Generally they bring 
in a different set of operator expectation values, but are fewer in number compared to those in the even chirality sum 
rules, at least to the leading order. It would be useful to consider both sets of sum rules for numerical evaluation. 

The behaviour of the absorptive parts in the space-like region in different sum rules is worth mentioning here. Such 
absorptive parts are a characteristic feature of thermal field theories and here it describes pion absorption by the 
current from the heat bath and its emission into it. For nonzero \q\, it is present in all the sum rules. However, as 
\q\ — > 0, it does or does not survive depending on the kinematics. In the sum rules considered here, it survives only 
in the sum rule (4.2) , obtained from the longitudinal amplitude defined in Ref. ||. 

We have evaluated the sum rules to order T 2 in the chiral limit. The sum rules are expected to be well saturated to 
this order and should be independent of the Borel parameter M 2 over a larger range. The subtraction eliminates the 
continuum contribution beyond about lGeV, it being practically independent of T for T < lbOMeV, say. There are 
thus no other significant source of contribution in the spectral representation besides those already included, namely, 
the 2ir -continuum and the p-resonance. 

As regards the operators, only those which are Lorentz scalars and appear already in the corresponding vacuum sum 
rules, may produce terms of order T 2 . We note, however, that while the T 2 contributions in the spectral representation 
arise only from 2tt and p states, scalar operators of all dimensions can in general contribute to this order, though, of 
course, the usual Borel suppression operates for operators of higher dimensions. 
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The other point to note is that the perturbative contribution, which is the coefficient of the unit operator, is 
eliminated by subtraction. So the higher dimension operators, which constituted only (power) corrections to the 
perturbative result in the vacuum sum rules, are now assuming the leading role. Thus the subtracted version of sum 
rules considered here brings out clearly the sensitivity of the thermal sum rules to any uncertainty in the operator 
expectation values. 

The present work is a continuation of our earlier work H . There our main concern was to include all the operators 
up to dimension four with their renormalisation effects in the vector- vector sum rules. But the omission of the scalar, 
four quark operator was not justified. For, while the dimension four operators are significant at higher temperatures 
near the critical point, at low temperatures it is only this dimension six operator which dominates, contributing to 
order T 2 in the chiral limit. 

The most extensive earlier work is that of Hatsuda et al 0. They include contributions up to order T 6 from 
the operators and plot the T-dependence of m p (T) and F p (T) ( and also parameters of other resonances) but do 
not obtain the coefficients of T 2 ,T 4 etc. (Incidentally their evaluation up to order T 6 is incomplete, because they 
include contributions up to order T 4 and T 6 only when they arc leading in some operators, but neglect non-leading 
contributions even to order T 4 in others, eg, (uu).) Eletsky et al (28) calculate the p meson mass to order T 4 , but do 
not take the mixing of dimension four, Lorentz non-scalar operators under renormalisation. 
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